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ABSTRACT 

We have compiled the largest sample of multiwavelength spectral energy distributions (SEDs) of 
Broad Absorption Line (BAL) quasars to date, from the radio to the X-ray. We present new Spitzer 
MIPS (24, 70, and 160/im) observations of 38 BAL quasars in addition to data from the literature 
and public archives. In general, the mid-infrared properties of BAL quasars are consistent with those 
of non-BAL quasars of comparable luminosity. In particular, the optical-to-mid-infrared luminosity 
ratios of the two populations are indistinguishable. We also measure or place upper limits on the 
contribution of star formation to the far-infrared power. Of 22 (57%) upper limits, seven quasars have 
sufficiently sensitive constraints to conclude that star formation likely contributes little (< 20%) to 
their far-infrared power. The 17 BAL quasars (45%) with detected excess far-infrared emission likely 
host hyperluminous starbursts with Lfir.SF = 1O^'^~^'*L0. Mid-infrared through X-ray composite 
BAL quasar SEDs are presented, incorporating all of the available photometry. Overall, we find 
no compelling evidence for inherent differences between the SEDs of BAL vs. non-BAL quasars of 
comparable luminosity. Therefore a "cocoon" picture of a typical BAL quasar outflow whereby the 
wind covers a large fraction of the sky is not supported by the mid-infrared SED comparison with 
normal quasars, and the disk- wind paradigm with a typical radio-quiet quasar hosting a BAL region 
remains viable. 

Subject headings: galaxies: active — quasars: absorption lines 



1. INTRODUCTION 

Recent theoretical investigations into quasar feedback 
as a potentially important cosmological force in galaxy 
formation (e.g., Scannapieco & Oh 2004; Granato et al. 
2004; Springel et al. 2005; Hopkins et al. 2005) have 
led to a expansion of interest in quasar outflows. En- 
ergetic outflows are most obviously manifested in Broad 
Absorption Line (BAL) quasars. This population, ap- 
proximately 20% (once selection effects are taken into 
account) of optically selected type 1 quasar samples (e.g., 
Hewett & Fohz 2003; Reichard et al. 2003a), is notable 
for broad, blueshifted absorption evident in common ul- 
traviolet resonance transitions such as C IV, Lya, and 
O VI. These P Cygni-type features arise because the ob- 
server is looking through an outflowing wind. Both the 
high velocities (up to 0.03 to 0.1c) and the relatively high 
ionization state of the gas indicate that the wind origi- 
nates within the inner parsec of the central engine, and 
is likely more or less co-spatial with the broad-line region 
gas (Murray et al. 1995; Proga et al. 2000). 

While the line of sight to an individual BAL quasar 
certainly cuts through an outflowing wind, the average 
covering fraction of the outflow remains difficult to con- 
strain. For example, it may be that BAL quasars have a 
wind that covers a large fraction of the sky, and therefore 
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only about 20% of quasars have BAL outflows. Alter- 
nately, most quasars may host BAL outflows, with the 
20% (corrected) population fraction resulting from the 
average BAL wind covering fraction. The latter interpre- 
tation is supported by evidence from spectropolarime- 
try that there are absorption-line-free lines of sight in 
BAL quasars (e.g.. Ogle et al. 1999), and the wind cov- 
ering fraction is inferred to be low from constraints on 
ultraviolet emission-line scattering (Hamann et al. 1993). 
Weymann et al. (1991) found an overall similarity of BAL 
and non-BAL quasar ultraviolet-optical continuum and 
emission-line properties. More subtly, the signiflcantly 
larger spectral comparison of Reichard et al. (2003b) re- 
vealed that BAL quasars are preferentially drawn from 
the non-BAL quasar population with intrinsically bluer 
continua, though they concluded that the parent popu- 
lations of optically selected BAL and non-BAL quasars 
appear to be the same. At submm wavelengths, where 
the thermal dust emission is generated on much larger 
scales than the accretion disk, the populations of BAL 
and non-BAL quasars have indistinguishable flux distri- 
butions, though the majority of both populations have 
only upper limits (Lewis et al. 2003; Willott et al. 2003; 
Priddey et al. 2007).^ These empirical results are consis- 
tent with the disk-wind paradigm whereby a radiatively 
driven, equatorial wind from the accretion disk gives rise 
to both the broad emission and broad absorption-line fea- 
tures (Murray et al. 1995). However, a near-equatorial 
line of sight has been called into question by radio obser- 
vations where typical orientation indicators (e.g., steep 
vs. flat radio spectra, and core vs. lobe-dominated flux 

® Priddey et al. (2007) also found the intriguing result that of 
their 16 BAL quasars, the submm-detected objects tended to have 
the largest C IV absorption EW. 



2 



Gallagher et al. 



ratios) suggest a range of inclination angles for radio- 
loud BAL quasars (Becker et al. 2000; Zhou et al. 2006; 
Brotherton et al. 2006). Furthermore, the identification 
of a few BAL quasars with ultra- luminous infrared galax- 
ies (ULIRGS; e.g., Canalize & Stockton 2001) has led to 
the interpretation that BAL quasars may be "cocooned" 
during an obscured evolutionary stage in the transition 
from galaxy merger to starbursting ULIRG to shrouded 
BAL quasar to naked quasar (e.g., Sanders et al. 1988; 
Gregg et al. 2002). 

The apparently contradictory results of these studies 
have led us to undertake a systematic survey of a well- 
defined sample of BAL quasars drawn from the Large 
Bright Quasar Survey (LBQS; Hewett et al. 1995; Hewett 
& Foltz 2003). The goal is to obtain multiwavelength 
coverage of a sizable number of BAL quasars for com- 
parison of the spectral energy distributions (SEDs) of 
BAL and non-BAL quasars in an effort to investigate 
the underlying relationship between the two populations. 
To date, there are two notable differences between BAL 
and non-BAL quasar SEDs: (1) BAL quasars typically 
show more evidence for dust reddening and extinction 
in their ultraviolet-optical spectra (Sprayberry & Foltz 
1992; Richards et al. 2002; Reichard et al. 2003b; Trump 
et al. 2006), and (2) BAL quasars are remarkably weak 
X-ray emitters (e.g., Green & Mathur 1996; Gallagher 
et al. 2006). However, both of these properties are con- 
sistent with orientation-dependent obscuration effects. 

In this present study, we focus on incorporating near 
and mid-infrared photometry into the picture, as longer 
wavelength data are expected to be less sensitive to ab- 
sorption and less influenced by orientation effects. In par- 
ticular, we will test the conclusion of Krolik & Voit (1998) 
that limb- darkening of the accretion disk emission should 
lead BAL quasars to have relatively weak ultraviolet- 
optical continua compared to non-BAL quasars if the 
BAL outflow is strongly equatorial. Goodrich (1997) 
reached a similar conclusion, though the proposed mech- 
anism was orientation-dependent attenuation (perhaps 
via electron scattering) based on polarization studies. 
The ultraviolet-optical emission in quasars is generated 
by the accretion flow near the black hole, while in- 
frared continuum light longward of ^ 1 /im is gener- 
ally attributed to a parsec-scale, cold, and dusty region 
that is heated by and reprocesses the higher energy, di- 
rect continuum. Therefore, the mid-infrared power in 
BAL quasars should provide an independent measure of 
the intrinsic accretion luminosity to test these hypothe- 
ses. 

Throughout we assume A-CDM cosmology with ilu = 
0.3, = 0.7, and Hq = 70kms-i Mpc-^ (Spergel et al. 
2003, 2007). 

2. SAMPLE SELECTION 

This sample was drawn from the Large Bright Quasar 
Survey (Hewett et al. 1995) sample of BAL quasars iden- 
tified by Hewett & Foltz (2003). Each object has a BAL 
probability of 1.0 and z > 1.4, the redshift at which 
the definitive C IV BAL is shifted into the wavelength 
regime accessible with ground-based spectroscopy. Of 
the 44 total BAL quasars meeting these criteria identified 
by Hewett & Fohz (2003), 38 were chosen for Spitzer 24, 
70, and 160/im MIPS observations; we present this MIPS 
sample in this paper. The MIPS targets were chosen 



to maximize the available multiwavelength coverage. In 
addition to the mid and far-infrared coverage, available 
data archives and literature were searched to augment 
the multiwavelength photometry as much as possible. 
The targets are optically bright with Bj^ magnitudes 
of 16.68-18.84, and they are drawn from a homogeneous, 
magnitude- limited quasar survey that has effective, well- 
defined, and objectively applied selection criteria (e.g., 
Hewett et al. 2001); six known BAL quasars with low- 
ionization (Mg II and/or Al III) BALs (hereafter, LoB- 
ALs) are included. The overall sample is luminous, with 
Mb of -26.2 to -28.3 (using the K-corrections for non- 
BAL quasars of Cristiani & Vio 1990). 

Our observed MIPS sample is listed in Table 1 includ- 
ing the redshift and BAL type. The classifications are 
HiBAL (no Al III or Mg II BAL detected), LoBAL (Mg II 
BAL detected), and Unknown (the electronic spectra do 
not cover Mg II). 

3. MULTIWAVELENGTH PHOTOMETRY 

The foundation of this multiwavelength SED project is 
complete mid and near-infrared photometry from Spitzer 
MIPS and the Two Micron AU Sky Survey (2MASS; 
Skrutskie et al. 2006) of 38 LBQS BAL quasars. These 
data have been supplemented by near-complete (35 of 
38) 0.5-8.0 keV X-ray data, multiband optical photom- 
etry for all, submm coverage for 16 objects, and Very 
Large Array (VLA) radio observations of 35 out of 38. 
In this section, the provenance of the SED data is briefly 
described. 

3.1. Spitzer MIPS Photometry 

The bulk of the MIPS data is the 36 BAL quasar obser- 
vations from the Cycle 1 Director's Discretionary Time 

for Gallagher's Spitzer Fellowship program. Two ad- 
ditional objects, 0059-2735 and 1331-0108, are added 
from Program #82 (PI Rieke) of Cycle 1 Guaranteed 
Time. 

The MIPS observation modes were one 3 s cycle at 
24/^m (48 s integration), one 3 s or 10 s cycle at 70/im 
(38 or 128 s integration), and four 10 s cycles at 160/im 
(84 s integration). The longer 70/im integration times 
were chosen for fainter sources, and the 160/im integra- 
tion times were chosen to reach the confusion limit. All 
BAL quasars were detected with high S /N at 24/im. We 
consider BAL quasars with S/N> 3 to be detected at 70 
and 160/Lim. The upper limits arc calculated as the flux 
in the source aperture -|- 3a or 3a, whichever is greater 
(following the convention of the MIPS instrument team), 
where a is the standard deviation in the mean back- 
ground. By this criterion, 13 BAL quasars are detected 
at 70/im, and three are detected at 160/im. (The two 
Guaranteed Time targets do not have 160/im coverage.) 
The Spitzer photometry is presented in Table 2. 

3.2. Near-infrared 

Thirty-two out of 38 BAL quasars are included in 
the 2MASS Point Source Catalog (PSC; Skrutskie et al. 
2006); of these, all had JHK detections except for 
00211-0213, 1133-h0214, 1314-^0116, and 2116-4439, 

^ This photographic blue magnitude is related to the more stan- 
dard Johnson blue magnitude: Bj = B — 0.28(-B — V) (Hewett 
et al. 1995). 
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which only have upper Umits in K . To obtain photom- 
etry (either upper limits or detections) for the remain- 
ing six (0051-0019, 1208-hl535, 1240+1607, 1243-^0121, 
1443+0141, and 2140-4552), we downloaded the best 
Atlas JHK images (recommended for photometry) from 
the 2MASS database.^ These images have all been flux- 
calibrated, combined, and resampled to a l"/pixcl scale; 
the photometric zeropoints are included in the fits im- 
age headers. We measured the signal (in data numbers) 
within a 2"-radius circular source cell centered on the 
LBQS optical position; the background was determined 
from an encircling annulus with inner and outer radii of 
8" and 10", respectively. The noise per pixel, a, is the 
standard deviation in the mean background. If the S/N 
within a 2" aperture was more than 2.5, we consider the 
source to be detected. Otherwise, the upper limit is set 
to the 3(7 background level within the source cell area. 

While a 2" source cell maximized the S/N, a signifi- 
cant fraction of the flux in the point spread function is 
distributed outside of this aperture. To determine the 
aperture correction (which varies from fleld to field and 
band to band), we performed aperture photometry on 
the 2MASS PSC catalog sources with 'AAA' photomet- 
ric quality flags and no galaxy contamination within 300" 
of the target that coincided with the Atlas image area. 
The aperture correction for each target in each band was 
then the mean difference between the catalog magnitudes 
and their corresponding 2"-aperture magnitudes. The 
number of point sources used in the aperture correction 
determination was between four and 16; aperture correc- 
tion values ranged from —0.42 to —0.57 magnitudes. The 
aperture correction was added to both the detections and 
upper limits to give the fluxes (converted to mJy from 
magnitudes) listed in Table 3. The standard deviation 
in the mean aperture correction has been folded into the 
photometric error. 

Following the procedure outlined above, we obtain 
37/35/32 detections in the J /H /K bands, respectively, 
for the 38 objects in the sample. 

3.3. X-ray 

Thirty-four of the 38 objects in the sample were 
observed with Chandra, and a detailed analysis of 
these data and comparison with ultraviolet spectro- 
scopic properties is presented in Gallagher et al. (2006). 
In addition, 2212—1759 was observed for 172 ks with 
XMM-Newton (Clavel et al. 2006). Using the 0.5- 
2.0 keV EPIC-pn count-rate upper hmit (< 9.25 x 
10~^ counts"^), and assuming a power-law spectral in- 
dex of ax = 0.0 (appropriate for absorbed sources), 
the modeling tool WebPIMMS'' gives /2keV < 4.12 x 

•|^Q-34 j,j.g g-1 {J2^^. 

In Table 3, the observed- frame 1 keV photometry is 
listed in units of 10~^ mJy for the 35 sources with X-ray 
data. The asymmetric flux errors result from Poisson 
X-ray counting statistics. 

3.4. Optical 

We have cross-correlated the LBQS Spitzer 
BAL quasar sample with the SDSS Data Release 5 
archive. Of the 38 MlPS-observed BAL quasars, 22 have 

* http;//irsa.ipac. caltech.edu/applications/2MASS/IM/ 
® http://heascirc.nasa.gov/Tools/w3pimms.html 



available SDSS photometry. We list the dereddened (for 
Galactic extinction) ugriz PSF magnitudes for these 
objects in Table 3 converted to units of mJy using the 
SDSS zeropoint of 3631 Jy. 

To obtain constraints on the optical colors of the re- 
maining quasars, we convolved the electronic spectra 
(from Korista et al. 1993 [preferred] or Hewett et al. 1995) 
with the SDSS gri filter functions for those quasars with- 
out SDSS coverage after first normalizing the available 
spectra using the LBQS Bj photometry and transmission 
curves (Appendix A; Maddox & Hewett 2006). These 
synthetic broad-band fiuxes are presented in Table 3 (in 
units of mJy), errors are propagated from the 0.15 mag 
Bj uncertainty (Hewett et al. 1995). The quoted uncer- 
tainty does not include additional errors from the con- 
volution. The number of bandpasses is limited to the 
wavelength coverage of the electronic spectra. 

To check on the absolute photometric accuracy of the 
synthetic fluxes, synthetic SDSS photometry was calcu- 
lated for those quasars with actual SDSS data using non- 
SDSS electronic spectra. In sum, the average offsets be- 
tween the synthetic and SDSS qri photometry (defined 
as (SDSS -''synthetic)/SDSS [mJy]) are within the 15% 
photometric uncertainty of the LBQS. The mean offsets 
and standard deviations of the mean offsets for ,9, r, and 
t are ± 27%, -5 ± 37%, and -12 ± 33%, respectively. 
The large standard deviations likely indicate that the 
synthetic photometric errors are underestimated, how- 
ever, there is no evidence for absolute, systematic offsets 
between the synthetic and SDSS photometry. 

3.5. Submillimeter 

Sixteen of the 38 BAL quasars have been observed at 
450 and 850/im with the Submillimetre Common-User 
Bolometer Array (SCUBA) instrument at the James 
Clerk Maxwell Telescope. The SCUBA data reduction 
and analysis is described in detail in Priddey et al. (2007). 
In this paper, we use S/N > 2.5 as the detection cri- 
terion; at 850/xm (450/xm), this yields five (two) detec- 
tions with a 1(7 limit of 1.5 mJy. A more conservative 
3(7 threshold reduces the number of 850/xm detections to 
three. For non-detections, the upper limits are given as 
the flux in the source aperture + 2a or 2a, whichever is 
greater. The available photometry is included in Table 2. 

3.6. Radio 

Radio coverage is publicly available for 35 of 38 
BAL quasars. These data are compiled from Stocke et al. 
(1992) (5 GHz), the Faint Images of the Radio Sky at 
Twenty cm (FIRST; White et al. 1997) (1.4 GHz), and 
the National Radio Astronomical Observatory VLA Sky 
Survey (NVSS; Condon et al. 1998) (1.4 GHz). Data are 
listed in Table 2. Upper limits for the FIRST survey are 
as given for the field in the catalog search; NVSS upper 
limits are 5 mJy (Condon et al. 1998). 

Only one of the quasars in our sample, 2211 — 1915, is 
known to be formally radio-loud as defined by a radio 
hmiinosity of > lO'^ergs^^ Hz^^. The three quasars 
without radio constraints are too far south for the VLA. 
Based on the radio-loud fraction of the LBQS determined 
from a cross-correlation with the FIRST survey (Hewett 
et al. 2001), we expect at most 12% (< 1) of these to 
be radio-loud, although a more likely value is 6% based 
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on the lower fraction of radio-detected BAL versus non- 
BAL quasars (Hewett & Foltz 2003). 

4. RESULTS AND DISCUSSION 

We present the available multiwavelength photometry 
for each BAL quasar in Tables 2 and 3. The SED of 

each quasar is plotted in Figures la-lf in rcst-framc units 
of log(z^/^) (ergs~^) vs. log{v) (Hz), where \og{liy) in- 
cludes the bandpass correction: log(/y) — log(l + z) + 
log(47r(i^^jjj) — 26 {f^ in mJy and luminosity distance, 
rfium, in cm). For reference, composite SEDs normal- 
ized to the 24/im MIPS data points are overlaid. The 
composite SEDs are the Elvis ct al. (1994) mean SED at 
log{v) < 13.57 (8^m) and the Richards et al. (2006) lumi- 
nous mean quasar SED for log(z^) = 12.50-17.00 (95/im 
to 0.4 keV). For the > 0.5 keV X-ray region, a spectral 
index of = —1.0 has been assumed (e.g., Gallagher 
et al. 2002; Page et al. 2005), with the normalization fixed 
by the best-fitting relation from Steffen et al. (2006, eq. 
Ic): log(Z2kev) = 0.721 log(Z25oo) +4.531. This accounts 
for the empirical observation that the most ultraviolet- 
luminous quasars emit relatively less power in X-rays. 
From a qualitative inspection of the SEDs, the most ob- 
vious discrepancy between the plotted c;omposites and 
the data is the known deficit at X-ray wavelengths when 
compared to the expected X-ray power. 

To investigate the nature of mid-infrared SEDs of 
BAL quasars, a comparison sample is necessary. Given 
the known SED changes of normal type 1 quasars with 
luminosity in the X-ray (e.g., Avni & Tananbaum 1986; 
Strateva et al. 2005; Steffen et al. 2006) and infrared 
regimes (Richards et al. 2006; Gallagher et al. 2007), 
identifying a comparison sample of objects with com- 
parable luminosity is clearly preferred. The largest sam- 
ple of quasar SEDs to date was presented by Richards 
et al. (2006) who incorporated the publicly available 
Spitzer MIPS-I-IRAG photometry as of the Data Re- 
lease 3 quasar catalog of Schneider et al. (2005). All 
259 objects were detected by Spitzer. We draw our 
comparison from the radio-quiet subset of this sam- 
ple matched in monochromatic 5000A luminosity to the 
LBQS BAL quasar sample (log(/^^^^^) = 31.21 32.29); 

a total of 40 objects. This luminous quasar sample is 
referred to as the SDSS comparison sample hereafter. 
While ideally the comparison sample would be drawn 
from the LBQS itself to mitigate possible selection ef- 
fects, sensitive multiwavelength data as required for the 
comparison do not exist elsewhere. Luminous quasars 
are quite rare, and so accumulating a sizable sample (of 
order tens) requires multiwavelength, large area surveys 
or explicit targeting. 

To compare directly the LBQS BAL quasar sample 
with the SDSS comparison sample, we compute three 
fiducial monochromatic luminosities at rest-frame 8/xm, 
5000A, and 2500A. The first is measured by normal- 
izing the Richards et al. (2006) luminous (log(iboi) > 
46.02 erg s~-^) quasar composite SED to the 24^m MIPS 
datapoint. The second two are calculated from the best- 
fitting normalization, Z,y,o, and spectral index, qiq, of a 
power-law model fit to the available rest-frame 1200- 
5000A photometry: li, = li,fiV°^° . The lower frequency 
bound for this ultraviolet-optical bandpass is chosen to 
avoid the optical spectral region of A > 5000A {v < 
•|^q14.78 JJ2;) where host-galaxy contamination can be sig- 



nificant (Vanden Berk et al. 2001); at A < 1200A(iy > 
-|^q15.4o jj^)^ the intervening Lya forest affects the spec- 
tra considerably. For the BAL quasars, the parameters, 
ao, log(/8Mm), log(/^^^^^), and log{l^^^^^), are listed in 
Table 1 along with log(Z2kcv) for those with X-ray data. 

For reference, two integrated infrared luminosities, 
Lir,Qso (1mm to 2/zm) and Lfir.QSO (1mm to 20/im), arc 
calculated using Igfim and the composite SEDs of Elvis 
et al. (1994) (for A > 40/ini) and Richards et al. 2006 
(for A < 40/im). Elvis et al. (1994) (which used IRAS 
data) has better coverage at longer wavelengths, though 
we point out that neither composite is empirically con- 
strained in the range from log(z^) = 10.2-12.5 (16 GHz 
to 95^m). Using the combined composite SED, the in- 
frared and far-infrared bolometric corrections from i^lsfj.m 
to iir.QSO and Lfir^Qso are 2.99 and 0.84, respectively. 
Given that the near and mid-infrared emission is believed 
to be dominated by the quasar, these integrated infrared 
luminosities (normalized at rest-frame ~ 8/im) may be 
considered an estimate of the quasar's contribution to 
the infrared SED. Any additional power in this regime 
could therefore be attributed to star-formation in the 
host galaxy. 

4.1. Ultraviolet- Optical Comparison 

Spectroscopic comparisons of BAL and non- 
BAL quasars include the seminal study by Weymann 
et al. (1991) and subsequent SDSS work by Reichard 
et al. (2003b). While Weymann et al. (1991) concluded 
that the underlying ultraviolet continuum and emission- 
line properties of BAL quasars are consistent with those 
of normal quasars from their sample of 42 BAL quasars, 
the larger (224) study of Reichard et al. (2003b) found 
that BAL quasars are more specifically drawn from a 
parent sample of intrinsically blue quasars. A further 
caveat is the statistical detection of excess dust extinc- 
tion and reddening in the BAL quasar population, most 
notably among LoBAL quasars (Sprayberry & Foltz 
1992; Trump et al. 2006). 

In this study, we evaluate the ultraviolet-optical con- 
tinua of BAL and non-BAL quasars in a complemen- 
tary fashion, using the rest-frame 1200--5000A photome- 
try. The inclusion of 2MASS data for the BAL quasars 
and the SDSS 2;-band photometry enables an extension 
of previous studies to the rest-frame optical. First, we 
compare the distribution of ao for both the LBQS BAL 
and SDSS comparison quasar samples; histograms of the 
two distributions are shown in Figure 2. 

To compare the two samples, we calculated the Stu- 
dent's T-statistic of ao, T = 0.30. This value indi- 
cates a probability of the null hypothesis that the two 
samples have the same mean of 77%. However, inspec- 
tion of Figure 2 reveals that this statistic is perhaps 
not capturing differences between the two samples: the 
BAL quasars seem to be slightly redder on average than 
normal quasars. The median values of -0.45 (SDSS) and 
0.63 (BAL quasar) reflect this apparent difference. To 
investigate if this is significant, we use the ASURV imple- 
mentation (La Valley et al. 1992) of two non-parametric 
statistical tests, the Wilcoxon and logrank, and obtain 
probabilities of the null hypothesis (i.e., that the BAL 
and SDSS comparison samples have the same distribu- 
tion of ao) of 0.04 and 0.08, respectively. These results 
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are not statistically significant evidence for a difference 
between the samples. 

To investigate further the distributions, we first ex- 
amine the effect of convolving a non-BAL quasar com- 
posite spectrum from the SDSS (Vanden Berk et al. 
2001) with the SDSS filter functions in the redshift 
range, z — 1.5-2.9, and then measuring cto- The strong 
Lya and C IV emission lines at the blue end of the 
1200-5000A bandpass give a bluer of -0.34 (aver- 
aged over the redshift range) than the value from fit- 
ting the continuum of the SDSS composite spectrum, 
3p^^=-0.44 (Vanden Berk et al. 2001). Using the SDSS 
BAL and LoBAL quasar composite spectra (Reichard 
et al. 2003a, b) in the same manner, we obtain ao=~0.67 
and -1.28, respectively. The BAL quasar composites give 
redder values of ao because of both dust reddening and 
extinction and the loss of flux on the; blue end from C IV 
and Si IV BALs. Both of these effects are stronger in 
the LoBAL quasars, and the typical extinction measured 
in Hi and LoBAL quasar spectra is E{B — V)=0.02 and 
0.08, respectively (Reichard et al. 2003b). This is con- 
sistent with the trend in the data where the median ao 
values are -0.59 and -0.90 for Hi and LoBAL quasars, 
respectively. We do not expect the composite values to 
represent our data precisely, because of the range of avail- 
able photometric data as well as the inherent dispersion 
in the ultraviolet-optical spectral properties of BAL and 
non-BAL quasars. 

The right-hand panel of Figure 2 illustrates an inter- 
esting artifact of optical quasar selection. In the plots 
of ao vs. / °, there is an apparent correlation be- 

5000A' 

twcx^ii tlic^ two whereby more optically luminous quasars 
are redder. This is a statistically significant trend in 
both samples, the values of Spearman's p of —0.498 and 
—0.452 indicate probabilities that the properties are not 
correlated of 0.0015 and 0.0034 for the BAL and SDSS 
comparison quasar samples, respectively. This is likely 
a selection effect, as the identical test of ccq vs. I 

' ° 2500A 

yields no evidence for a significant correlation with prob- 
abilities of 0.140 and 0.107 for the respective samples. 
The monochromatic luminosity at 2500A is closer to the 
selection bandpasses of Bj (for the LBQS) and i (for the 
SDSS) in this redshift range. Of two quasars with similar 
fluxes in the selection bandpass, the redder one will be 
more luminous at longer wavelengths. 

Overall, therefore, we find no compelling evidence for 
differences in the ultraviolet-optical continua in BAL 
and non-BAL quasars beyond the known tendency of 
BAL quasars to be mildly dust-reddened, in agreement 
with previous optical spectroscopic studies. As the opti- 
cal quasar selection criteria of the SDSS and the LBQS 
are different, there is some concern that selection effects 
could be important. Briefly, the SDSS selection for radio- 
quiet quasars is based on an i-band flux limit and broad- 
band photometric color selection, whereas the LBQS has 
a Bj flux limit and a complicated selection algorithm 
using objective prism spectra. LBQS quasar candidates 
were identifled via blue spectra and/or from features such 
as emission lines and continuum breaks (§3.1.4; Hewett 
et al. 1995). Though the SDSS should find aU LBQS 
quasars, the converse is not necessarily true. (An ex- 
plicit comparison with the FIRST Bright Quasar Sur- 
vey [White et al. 1997] determined an incompleteness 



of 13±4% for the LBQS [Hewett et al. 2001]). There- 
fore, it is reassuring that only four of the 40 SDSS com- 
parison quasars have values of ao redder than those of 
the BAL quasars. In subsequent analyses we examine 
whether these four objects affect our results. 

4.2. Mid-Infrared Comparison 

Prom visual inspection of the mid-infrared SEDs pre- 
sented in Figure 1, the composite quasar SED normal- 
ized at the 24/im MIPS data point appears to predict 
quite well the ultraviolet-optical photometry for most 
BAL quasars. As there is an inherent dispersion of 
quasar SEDs that is not captured in the mean compos- 
ite (Elvis et al. 1994; Richards et al. 2006), we quan- 
tify the relationship between the mid-infrared and op- 
tical by comparing /g/tm and l^^^^^J^ for the BAL and 
SDSS comparison samples. Specifically, we measure lin- 
ear least-squares fits for log(Z8^ji„) vs. log(Z^^^^^^^^ ^) and 

log(/ "/Zsum) VS. log(Z '). As can be seen in Fig- 

°^ 5000A' ' °^ 5000A' ° 

ure 3, both samples occupy the same region of parameter 

space, and the normalizations and slopes arc consistent 
within la. Note that by using / ° rather than I 

° 5000A 2500A' 

we mitigate the effects of optical extinction; at 5000A, 
the monochromatic luminosity is only reduced by ~ 7% 
for E{B -V) = 0.08. 

In addition to consistent fits with the data, we also 
measure the Student's T-statistic of log(/ ° /hum) for 

'^^ 5000A' ' 

both samples; T = —1.39. This value indicates a 17% 
probability that the null hypothesis (that the two sam- 
ples have the same mean) is correct. Inspection of his- 
tograms of the two distributions and similar medians 
(-1.09/-1.14forSDSS/BAL) confirm this conclusion. Fi- 
nally, the Wilcoxon and logrank tests give null hypoth- 
esis probabilities of 16% and 17%, respectively. If the 
four SDSS comparison quasars with values of ao redder 
than the BAL quasars are excluded, the Student's T-test 
significance increases to 73% that the samples have the 
same mean, and the medians become almost identical: 
-1.13 versus -1.14. Therefore, we find no evidence of a 
discrepancy in mid-infrared relative to optical power for 
BAL and non-BAL quasars of comparable luminosity. 

4.3. Far-Infrared Properties 

Prom inspection of Figure 1, the far- infrared (A < 
20/xm) properties of the BAL quasars arc not homoge- 
neous. Specifically, if one focuses on the quasars with 
70, 160^m and/or submm detections (rather than just 
upper limits), there is a notable dispersion in relative 
and absolute far- infrared power. This is in contrast 
to the mid-infrared through optical regimes, where the 
BAL quasar SEDs much more closely follow the compos- 
ite SED. Though most of the 70 and leO/xm MIPS and 
SCUBA data points are upper limits, several are sensi- 
tive enough to put useful limits on far-infrared emission 
in excess of a composite quasar SED. To estimate the 
possible contribution of star formation, we subtract the 
composite; quasar SED normalized to the 24//m MIPS 
point from the far-infrared photometry. An illustrative 
starburst SED from Efstathiou et al. (2000) is then nor- 
malized to the excess emission and integrated from 1mm 
to 20/^m to obtain igr^sF- The starburst model, char- 
acterized by a burst age of 16 Myr and a visual opti- 
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cal depth of 150, was generated for the prototype far- 
infrared/submm- hyper luminous z > 4 radio-quiet quasar 
BR1202— 0725. This object is sampled densely enough in 
the far-infrared through submm region that, unlike the 
majority of high redshift, submm-luminous quasars, well- 
constrained model fits to its photometric data can be ob- 
tained (Priddey & McMahon 2001). Ideally, one would 
prefer to generate individual fits to the far-infrared data, 
however, the data quality available here does not warrant 
it. 

The ifir.SF values are listed in Table 1, and the star- 
burst models are overplotted in Figure 1. Upper limits 
are given when there are no far-infrared detections or the 
data are consistent with a quasar SED. In the latter case, 
i'fir.SF was set to 0.2Lfir,QSO to take into account reason- 
able photometric uncertainties - the best-fitting values 
were significantly lower. 

In Figure 4, the estimated contributions from star 
formation are plotted against the quasar far-infrared 
power. Of 38 BAL quasars in the sample, 22 (58%) 
have only upper limits to ifir,SF, of these, seven (18%) 
have far-infrared constraints sufficient to determine that 
the quasar likely dominates their far-infrared emission. 
Based on local examples of a few LoBAL quasars 
identified with ULIRGS, an anecdotal identification of 
LoBAL quasars with merger remnants and starburst ac- 
tivity had been suggested (e.g. Low et al. 1988; Canal- 
izo & Stockton 2001). Though the two most Imni- 
nous starbursts are the LoBAL quasars 1331—0108 and 
2358-1-0216, the other four are not obviously different 
from the HiBAL quasar population. However, the pre- 
ponderance of upper limits at far-infrared wavelengths 
for the majority of Lo and HiBAL objects hampers our 
ability to evaluate this topic further. 

Of the 17 (45%) BAL quasars with Lfir,sF detec- 
tions, their Lfir.sF values range from 1.4-13ifir,Qso or 
IO^^^^^Lq, qualifying them as hyperluminous starbursts. 
Using the standard conversion^" to star formation rates, 
SFR « lO-iOLfir,sF/i-0 (M0 yr-i), gives SFR = lO^-^ 
Mq yr~^. Given the upper limits, even the objects with- 
out a detected far-infrared excess could host substantial 
starbursts. However, the uncertain origin of far-infrared 
emission in quasars, from cold dust heated by star forma- 
tion and/or accretion power (see Rowan- Robinson 2000 
and Haas et al. 2003 for opposite views on this matter), 
means that a far-infrared 'excess' could be accretion pow- 
ered. 

4.4. A Mid-Infrared-to-X-ray BAL Quasar Composite 

SED 

To synthesize the results presented in §4 thus far, we 
construct three composite SEDs, one for the SDSS com- 
parison sample, a second for the entire BAL quasar sam- 
ple of 38, and a third for the 25 known HiBAL quasars. 
The composite from the SDSS comparison sample has 
been made with photometric data (including IRAC) from 
the 36 quasars whose values of Qq are greater (bluer) than 
— 1.45, the minimum for the BAL quasars. With only 
six LoBAL quasars, there were not enough data to con- 
struct a composite SED for this type. The photometry 
from 24/im through 1 keV is used; there are insufficient 

The uncertainty in this relation is dominated primarily by the 
assumed initial mass function (cf. Priddey & McMahon 2001). 



detections at longer wavelengths to incorporate the far- 
infrared and submm. 

As a first step, we normalize the SEDs for each ob- 
ject to (at log(^')=14.78); this is necessary to use 

the redshift spread (z = 1.5 2.9) to fill out the multi- 
wavelength wavelength coverage. Next, we use a vari- 
able frequency window sized to contain at least seven 
photometric data points (though a few bins at the high 
frequency ends of available photometric coverage have 
fewer data points). Within each bin, the median nor- 
malized luminosity is identified. In Figure 5, we plot the 
normalized data for all BAL quasars, as well as our SDSS 
and BAL quasar composites. The Richards et al. (2006) 
non-BAL quasar composite is overplotted for reference. 
The X-ray regime of the SDSS composite has been ap- 
propriately normalized to match the median log(^ ° ) 

2500 A. 

(31.42 erg s~^) BAL quasar luminosity as described in 
the beginning of §4. 

From the BAL quasar composite, the mid-infrared 
emission shows no evidence for being brighter than ex- 
pected based on the 5000A luminosity in comparison 
with the SDSS composite. In the optical-ultraviolet re- 
gion, the structure of the composites follows the bright 
emission lines, including Ha, H/?, Mg II, C IV, and 
Lya. The reddened continuum in comparison with the 
Richards et al. (2006) composite is also evident, as is 
very strong X-ray absorption. While the BAL and Hi- 
BAL quasar composites are not noticeably different, the 
LoBAL quasar data cluster in the lower panel of Figure 5. 
In particular, the ultraviolet-optical photometry indi- 
cates stronger dust reddening and extinction and the ex- 
treme X-ray weakness is pronounced. The mid-infrared 
behavior bears further investigation, as some of the 24/im 
data appear to lie at the faint end of the distribution. 
The composite SDSS comparison and BAL quasar SEDs 
are listed in Table 4. 

5. SUMMARY AND CONCLUSIONS 

Using existing archives and literature as well as new 
Spitzer MIPS observations, we present 38 BAL quasar 
spectral energy distributions, the largest sample of ra- 
dio through X-ray data of BAL quasars to date. We 
compare this sample to a radio-quiet SDSS quasar sam- 
ple matched in 5000A luminosity whose data were first 
presented by Richards et al. (2006). In agreement with 
the conclusions of Weymann et al. (1991) and Reichard 
et al. (2003b), we find that luminous, optically selected 
BAL and non-BAL quasars have consistent ultraviolet- 
optical continua once dust reddening is taken into ac- 
count. In addition, the relative mid-infrared and optical 
power in both populations are indistinguishable. As al- 
ready known, the most notable difference between BAL 
and non-BAL quasar SEDs is in the X-ray regime, where 
BAL quasars are weak X-ray emitters as a result of intrin- 
sic absorption (e.g., Gallagher et al. 2006). All of these 
characteristics are evident in the composite BAL quasar 
and HiBAL quasar SEDs presented in Figure 5 and Ta- 
ble 4. 

Quasar mid-infrared emission is believed to arise on 
spatially large (> 1 pc) scales by reprocessing of the 
ultraviolet-optical continuum emission by a cold, dusty 
medium. In the picture where the covering fraction of a 
dusty BAL wind is large, one might expect BAL quasars 
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to be relatively brighter in the mid-infrared than non- 
BAL quasars with little or no wind because of a larger 
emitting volume of dust. For reference, if all (i.e., assum- 
ing a 47r covering fraction for the dust) of the 1200-5000A 
continuum luminosity absorbed by dust with the typical 
extinction found for Hi and LoBAL quasars {E{B — V) = 
0.023 and 0.077, respectively; Reichard et al. 2003b) 
were reemitted in the mid-infrared. Hi (Lo) BAL quasars 
should be 18% (47%) brighter than typical quasars. This 
picture is not upheld by our data. 

An understanding that BAL quasars seen along equa- 
torial lines of sight should have preferentially fainter and 
redder ultraviolet-optical continua than normal quasars, 
by the mechanism of accretion-disk limb darkening (Kro- 
lik & Voit 1998) or orientation-dependent, gray atten- 
uation (Goodrich 1997), leads to a similar expectation 
of mid-infrared bright BAL quasars. Again, our results 
do not support this. However, the characterization of 
disk-wind models as generating 'equatorial' outflows is 
perhaps taken too literally; the lines of sight closest 
to the accretion disk are likely completely blocked by 
the hot dust that obscures the inner accretion region 
in type 2 (narrow-line) quasars and generates the mid- 
infrared emission. Furthermore, radiatively driven winds 
are radial with the inclination angle of the outflow de- 
pendent on the launching radius: winds launched from 
closer to the black hole have a smaller inclination angle 
with respect to the disk normal (Proga et al. 2000; Ev- 
erett 2005). Therefore, even in the disk- wind paradigm 
where limb-darkening and orientation-dependent attenu- 
ation are operating, it is unclear that the inclination an- 
gles of BAL and non-BAL quasars would be sufficiently 
different to detect a mid-infrared excess with the data in 
hand. 

Based on the comparison of mid-infrared and optical 
properties of BAL vs. non-BAL quasars, well-designed, 
wide-angle, mid-infrared quasar surveys offer the most 
promise for capturing the 'true' BAL quasar fraction of 
type 1 quasars. Near-infrared surveys sampling the rest- 
frame optical may also prove fruitful in the future, but 
are currently not nearly sensitive enough. Optically se- 
lected surveys will inevitably be somewhat incomplete 
because the ultraviolet BALs along with any dust extinc- 
tion conspire to make BAL quasars fainter than normal 
quasars of comparable intrinsic luminosity. Clearly, soft- 
X-ray surveys will suffer from similar (though much more 
extreme) problems because of the severe effects of intrin- 
sic absorption. At higher X-ray energies, absorption be- 
comes much less of a factor. However, the potential of 
hard (> 10 keV) X-ray surveys to obtain similar com- 
pleteness to Spitzer will not be realized for many years 
- 2; ~ 2 quasars are not within the flux limits of current 



Avni, Y. & Tananbaum, H. 1986, ApJ, 305, 83 
Barvainis, R. 1987, ApJ, 320, 537 

Becker, R. H., White, R. L., Gregg, M. D., Brothcrton, M. S., 
Laurent-Muchlciseii, S. A., & Arav, N. 2000, ApJ, 538, 72 

Brotherton, M. S., dc Brcuck, C, & Schacfcr, J. J. 2006, MNRAS, 
L90 

Canalizo, G. & Stockton, A. 2001, ApJ, 555, 719 

Glavel, J., Schartel, N., & Tomas, L. 2006, A&A, 439 

Gondon, J. J., Gotten, W. D., Greisen, E. W., Yin, Q. F., Perley, 

R. A., Taylor, G. B., & Broderick, J. J. 1998, AJ, 115, 1693 
Gristiani, S. & Vio, R. 1990, A&A, 227, 385 



facilities. In any case, there is evidence for Compton- 
thick (A''h> 1.5 X 10^^ cm^^) X-ray absorption in some 
BAL quasars that will affect even > 10 keV X-rays (e.g., 
Gallagher et al. 2006). 

Given the dearth of far-infrared coverage (particularly 
detections) of non-BAL quasars of similar redshifts and 
luminosities to this sample, it is unclear how the far- 
infrared properties of BAL quasars would fit within a 
matched type 1 population. At lower redshifts and lu- 
minosities, the contribution of star formation to the far- 
infrared power of quasars is quite variable, and many 
(non-BAL) quasars show evidence for quite luminous 
starbursts (Fritz et al. 2006; Schweitzer et al. 2006). Even 
though Spitzer is significantly more sensitive than IRAS , 
the confusion limit for MIPS is typically brighter than 
the expected 70 and 160/xm fluxes for luminous, 2; ~ 2 
quasars. This is unfortunate, given that the far-infrared 
is one spectral region likely to show conclusive evidence 
of a star-forming host galaxy. We draw attention, how- 
ever, to the value of submm data to set important con- 
straints on excess far-infrared emission; six of the seven 
BAL quasars with the tightest upper limits to ifir.sF have 
sensitive SCUBA data. Thus the advent of the next gen- 
eration submm and far-infrared detectors will shed new 
light on this issue. 

Though future studies with IRAC and IRS on Spitzer 
may reveal subtle differences in BAL vs. non-BAL quasar 
1 8/im continuum properties, the gross, intrinsic proper- 
ties of the SEDs of the two populations are fully consis- 
tent. Two regimes where they do differ, in the tendency 
to exhibit optical dust reddening and strong X-ray ab- 
sorption, can both be accounted for with orientation- 
dependent absorption in a wind; the ultraviolet-optical 
and X-ray continua are not intrinsically different. To 
date therefore, the disk-wind paradigm where all quasars 
have a BAL region, and the detection of such is merely 
an accident of orientation, has not been overturned. 
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46.53 


< 45.83 


1239+0955 


2.013 


Hi 


-0.59 ± 0.22 


26.20 


31.51 


31.69 


32.83 


46.88 


46.33 


< 45.63 


1240+1607 


2.360 


Hi 


-0.72 + 0.18 


25.77 


31.34 


31.55 


32.57 


46.62 


46.07 


< 47.18 


1243+0121 


2.796 


? 


-0.66 ± 0.08 


26.10 


31.71 


31.91 


33.31 


47.36 


46.81 


< 46.11 


1314+0116 


2.686 


7 


-0.87 + 0.05 


25.79 


31.65 


31.91 


32.67 


46.72 


46.17 


< 47.10 


1331-0108 


1.881 


Lo 


-1.45 ± 0.26 


25.57 


31.61 


32.05 


33.03 


47.08 


46.53 


47.43 


1442-0011 


2.226 


Hi 


-0.39 + 0.22 


< 25.78 


31.52 


31.64 


33.05 


47.10 


46.54 


< 47.34 


1443+0141 


2.451 


7 


-0.52 ± 0.20 


25.93 


31.48 


31.64 


32.99 


47.03 


46.48 


< 45.78 


2111-4335 


1.708 


Hi 


-0.12 + 0.02 


26.14 


31.83 


31.87 


33.00 


47.05 


46.50 


< 45.80 


2116-4439 


1.480 


Hi 


-0.32 ± 0.13 


25.00 


31.34 


31.44 


32.63 


46.68 


46.13 


46.37 


2140-4552 


1.688 


Hi 


-0.13 ± 0.04 


26.05 


31.17 


31.21 


32.60 


46.64 


46.09 


< 46.68 


2154-2005 


2.035 


Hi 


-0.51 + 0.09 


26.10 


31.45 


31.61 


32.97 


47.02 


46.47 


< 45.77 


2201-1834 


1.814 


Hi 


-0.59 + 0.60 


< 25.71 


31.63 


31.80 


32.65 


46.70 


46.15 


46.74 


2211-1915 


1.952 


Hi 


-0.71 ± 0.10 


26.82 


31.53 


31.74 


32.62 


46.67 


46.12 


46.56 


2212-1759 


2.217 


Hi 


-0.63 + 0.07 


< 24.68 


31.64 


31.83 


32.83 


46.88 


46.33 


47.18 


2350-0045A 


1.624 


Lo 


-0.63 + 0.11 


< 25.21 


31.07 


31.26 


32.54 


46.59 


46.04 


< 46.70 


2358+0216 


1.872 


Lo 


-1.37 + 0.49 


< 25.65 


31.47 


31.88 


32.78 


46.83 


46.28 


47.39 



^'Redshifts from Hcwctt & Foltz (2003). 

''Key: Hi = ultraviolet spectra show high-ioiiizatioii BALs only; Lo = low-ionization (Al III and/or Mg II) BALs; ? = BAL type unknown 
because of redshift (Lamy & Hutscmckcrs 2004; Gallagher ct al. 2006; Clavcl ct al. 2006). 
'^Spectral index for 1^ oc u"° from a fit to the rest-frame 1200-5000A photometry. 

"^Monochromatic 2 kcV luminosity converted from the Chandra data originally presented in Gallagher ct al. (2006). 

''The 2500 and 5000A monochromatic luminosities are derived from the power-law model fit to the rest-frame 1200-5000A photometry. 
*The 8^m monochromatic luminosity is calculated from the Richards et al. (2006) composite luminous quasar SED normalized to the 
24/im MIPS data point. 

^Spectrally integrated infrared (1mm to 2^m) and far-infrared (1mm to 20/im) quasar luminosities are calculated from i^is^im as described 
in §4. 

'^Spectrally integrated far-infared (1mm to 20^m) luminosity from star formation calculated by normalizing a starburst model to the 
far-infared emission in excess of a composite quasar SED (see §4.3). 



Table 2. Infrared, Submillimeter, and Radio Photometry'' 



Name Spitzer MIPS SCUBAt^ Radio 





24//ni 


7n/(m 
I U/AIli 








5 GHz*^ 


1 4 GHz*^ 


0004+0147 


7.05 + 0.13 


< 23.34 


< 143.88 






< 0.33 


< 5.00 


0009+0219 


1.80 + 0.16 


< 29.44 


< 158.84 








< 5.00 


0019+0107 


8.74 ± 0.18 


< 26.15 


< 220.99 


50.0 + 16.0 


8.2 + 2.3 


< 0.27 


< 0.92 


0021-0213 


5.57 + 0.15 


< 33.23 


< 203.47 


< 17.4 


5.3 + 1.1 


< 0.30 


< 1.00 


0025-0151 


3.58 ± 0.14 


< 35.68 


< 159.29 


< 32.0 


3.5 + 1.4 


< 0.36 


< 1.02 


0029+0017 


6.14 + 0.17 


28.18 + 7.73 


< 99.77 


< 41.0 


5.2 + 2.0 


< 0.27 


< 0.95 


0051-0019 


8.00 + 0.15 


< 25.30 


< 129.65 








< 0.92 


0054+0200 


5.12 ± 0.15 


< 20.67 


< 84.68 








< 5.00 


0059-2735 


11.27 + 0.29 


34.07 + 7.24 








< 0.36 


< 5.00 


0106-0113 


5.83 ± 0.15 


< 29.81 


< 497.86 








< 0.96 


0109-0128 


6.88 ± 0.17 


< 35.20 


< 463.49 








< 0.97 


1029-0125 


4.62 ± 0.18 


41.27 + 7.98 


< 171.81 


< 17.2 


< 3.4 


< 0.27 


< 1.14 


1133+0214 


8.78 ± 0.18 


< 27.83 


< 139.10 








< 1.01 


1203+1530 


4.10 + 0.16 


< 54.21 


< 88.74 








< 0.95 


1205+1436 


4.93 ± 0.16 


53.69 ± 10.54 


< 137.77 








< 0.96 


1208+1535 


4.78 ± 0.17 


43.79 ± 8.08 


< 97.89 


< 34.0 


< 5.0 


< 0.33 


< 0.95 


1212+1445 


3.80 ± 0.17 


22.28 + 6.65 


< 106.44 






0.46 + 0.08 


< 0.95 


1216+1103 


8.78 ± 0.17 


< 40.16 


< 100.47 






< 0.36 


< 0.95 


1224+1349 


7.35 + 0.15 


40.79 + 8.09 


< 117.84 








< 0.97 


1231+1320 


2.82 + 0.21 


30.13 + 8.31 


< 200.70 


< 16.3 


< 2.8 


< 0.24 


< 1.97 


1235+0857 


10.49 ± 0.14 


< 47.89 


< 89.45 


< 13.8 


< 3.4 


0.51 + 0.07 


< 1.02 


1235+1453 


5.63 + 0.15 


< 26.15 


< 76.00 


< 40.0 


< 4.7 


< 0.39 


< 0.91 


1239+0955 


6.90 ± 0.14 


< 40.33 


< 178.90 


< 10.8 


< 2.0 




< 0.98 


1240+1607 


2.65 + 0.15 


< 18.20 


< 120.86 






< 0.30 


< 0.98 


1243+0121 


9.89 + 0.18 


< 26.11 


< 193.57 


< 14.4 


< 3.8 


< 0.30 


< 1.05 


1314+0116 


2.50 + 0.14 


< 35.06 


< 87.90 






< 0.27 


< 0.90 


1331-0108 


12.89 + 0.37 


76.59 ± 7.97 








1.54 + 0.19 


2.99 + 0.14 


1442-0011 


9.04 + 0.15 


< 50.13 


< 174.74 






0.22 + 0.06 


< 0.95 


1443+0141 


6.32 ± 0.13 


< 23.75 


120.00 ± 36.73 


33.8 + 8.5 


5.2 + 1.2 


< 0.27 


< 1.00 


2111-4335 


14.81 ± 0.12 


< 26.51 


117.20 ± 34.82 










2116-4439 


8.80 + 0.12 


29.58 + 5.17 


< 90.92 










2140-4552 


6.01 ± 0.14 


< 20.70 


< 198.78 










2154-2005 


9.34 + 0.15 


< 26.93 


< 102.33 


< 32.0 


< 5.2 




< 5.00 


2201-1834 


5.78 ± 0.17 


30.73 + 9.95 


< 206.85 


< 26.0 


< 6.0 


< 0.33 


< 5.00 


2211-1915 


4.56 ± 0.14 


< 25.49 


69.62 ± 17.08 


< 38.0 


< 4.1 




64.00 + 2.00 


2212-1759 


5.59 ± 0.14 


31.33 + 8.93 


< 98.11 






< 0.33 


< 5.00 


2350-0045A 


5.81 + 0.14 


< 31.38 


< 125.70 






< 0.36 


< 1.00 


2358+0216 


7.31 + 0.21 


79.41 ± 10.49 


< 239.25 








< 5.00 



''AH flux densities are in units of mjy. 

''SCUBA data were originaUy presented in Priddey et al. (2007). See §3.5 for details on detection criteria and upper limits. 
<=AU the 5 GHz data are from Stocke et al. (1992). 

<*Data are compiled from the FIRST (White et al. 1997) and NVSS (Condon et al. 1998) catalogs. 



Table 3. X-ray, Optical and Near-Infrared Photometry*' 



(LBQS B) 


Chandra^ 
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K 
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it 
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± 
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it 


0.055 
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it 0.063 
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it 
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it 
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373 
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it 
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294 
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± 
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it 
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it 0.067 
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it 
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it 
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< 
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± 
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± 
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it 
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153 
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003 
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± 
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it 
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it 
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it 
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2 
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± 
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± 





002 
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it 


0.006 
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it 


0.062 





315 
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1 
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030 
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037 
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it 
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u 
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U 
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1 
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u 
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it 


n 
U 


nQQ 


n 

u 
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it 
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n 
U 
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'-''-'^ — 0.214 





200 


± 





028 
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1 
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004 
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1 
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233 
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003 





292 


1 





006 





297 


1 

± 


0.006 





355 


1 

± 





039 





320 


it 


0.063 





497 


it 0.070 


1029-0125 


1 


431+°-5^5 





181 


± 





025 















168 


± 





023 





219 


± 





030 





293 


± 





041 













379 


it 





042 





434 


± 


0.051 





584 


± 0.103 


1133-1-0214 


1 


540+°-^33 





189 


± 





026 





184 


lb 





004 





249 


it 





003 





281 


it 





004 





338 


it 





006 





327 


lb 


0.008 





296 


it 





045 





466 


it 


0.075 




< 


0.196 


1203-1-1530 





^f^o+0.164 
0.098 





141 


± 





020 





090 


± 





002 





126 


it 





002 





166 


it 





003 





216 


it 





004 





224 


lb 


0.005 





307 


it 





044 





382 


it 


0.064 





364 


it 0.057 


1205-1-1436 


IC 


400+l?i2 





189 


± 





026 





165 


± 





004 





198 


it 





004 





212 


it 





004 





269 


it 





004 





247 


lb 


0.006 





334 


it 





038 





256 


it 


0.064 





370 


it 0.055 


1208-1-1535 





762+°-^^2 

' "-^—0.297 





287 


± 





040 





137 


± 





004 





178 


it 





003 





220 


it 
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253 


it 





005 





283 


lb 


0.006 





189 
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046 




< 
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±0.061 


1212-1-1445 




< 0.463 





303 
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042 





188 


± 





005 
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455 


it 





007 





445 


lb 


0.008 





474 
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it 


0.069 
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it 0.062 


1216-1-1103 
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008 





274 


it 
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378 


it 
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lb 


0.007 
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it 
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it 
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433 


it 0.098 
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± 





031 





224 
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357 


it 
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378 


lb 
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462 
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048 
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it 
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''AH flux densities are in units of mjy, except for X-ray fluxes which are in units of 10 *' mjy. 

^Chandra X-ray data were all originally presented by Gallagher et al. (2006) with the exception of 2212 — 1759 observed with XMM-Newton and presented by Clavel et al. (2006). 
"^Photometry from the original LBQS survey. Typical survey errors of 0.15 mag have been converted to mjy (Hewett et al. 1995). 

"^SDSS data are dereddened Data Release 5 PSF magnitudes (York et al. 2000; http://www.sdss.org/dr5). Data with incomplete coverage are synthetic flux densities calculated by 
convolving the SDSS gri filter functions with electronic spectra (see §3.4 for further details). 
''Photometry from the 2MASS catalog (Skrutskie et al. 2006) or aperture photometry as described in §3.2. 



Table 4. Composite SEDs 



SDSS Comparison Quasars All BAL Quasars HiBAL Quasars 





lu/l : 

' 6000 A 


Number'^ 




l^/l ' 

' 6000 A 


Number'^ 




' 6000A 


Number'^ 


13.525 


1.261 


8 


13.509 


1.210 


8 


13.512 


1.210 


8 


13.572 


1.126 


7 


13.532 


1.225 


7 


13.548 


1.249 


7 


13.624 


1.040 


10 


13.565 


1.136 


9 


13.591 


1.178 


8 


13.710 


1.046 


8 


13.607 


1.114 


8 


13.619 


1.123 


2 


13.809 


0.944 


7 


13.664 


1.173 


6 


14.558 


-0.062 


8 


14.012 


0.856 


8 


14.555 


0.004 


8 


14.594 


0.042 


7 


14.060 


0.636 


7 


14.578 


-0.024 


7 


14.640 


0.053 


10 


14.110 


0.598 


10 


14.611 


0.011 


9 


14.674 


0.024 


7 


14.149 


0.652 


9 


14.646 


0.069 


7 


14.705 


-0.038 


8 


14.189 


0.502 


9 


14.674 


0.054 


11 


14.748 


-0.034 


7 


14.216 


0.384 


7 


14.702 


-0.038 


12 


14.787 


-0.027 


8 


14.240 


0.399 


8 


14.731 


-0.044 


9 


14.819 


-0.038 


7 


14.260 


0.469 


7 


14.770 


0.032 


8 


14.861 


-0.071 


7 


14.282 


0.107 


7 


14.799 


-0.061 


7 


14.900 


-0.077 


7 


14.314 


0.168 


9 


14.817 


-0.005 


9 


14.952 


-0.116 


7 


14.345 


0.200 


11 


14.838 


-0.079 


7 


15.017 


-0.086 


9 


14.364 


0.141 


9 


14.861 


-0.070 


7 


15.051 


-0.066 


7 


14.392 


0.001 


7 


14.899 


-0.077 


7 


15.086 


-0.198 


7 


14.419 


0.011 


8 


14.943 
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15.109 
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14.453 
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11 


14.981 
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7 


14.501 
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7 


15.024 
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8 


15.168 
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14.546 
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15.045 
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7 


15.199 
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7 


14.629 
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15.071 
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14.743 
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-0.080 


7 
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-0.094 


9 
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-0.257 


8 


15.315 
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15.037 


-0.071 


7 


15.199 


-0.245 


8 


15.342 


-0.354 


8 


15.056 


-0.086 


7 


15.225 


-0.276 


17 


15.409 


-0.533 


7 


15.080 


-0.169 


9 


15.240 


-0.136 
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17.797 


-5.900 


7 


15.113 


-0.162 
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15.256 


-0.356 


12 


17.835 


-5.469 


7 


15.135 


-0.183 


10 


15.274 


-0.367 


10 


17.875 


-5.666 


7 


15.154 


-0.173 


7 


15.290 


-0.414 


7 


17.902 


-5.727 


3 


15.187 


-0.245 


10 


15.315 


-0.331 


11 








15.218 


-0.173 


12 


15.338 


-0.357 


7 








15.235 


-0.139 


7 


15.350 


-0.430 


7 








15.268 


-0.259 


9 


15.374 


-0.359 


8 








15.295 


-0.254 


7 


15.395 


-0.446 


7 








15.315 


-0.226 


8 


15.463 


-0.777 


7 








15.337 


-0.258 


7 


17.795 


-6.102 


7 








15.361 


-0.222 


8 


17.818 


-5.651 


7 








15.391 


-0.276 


9 


17.854 


-5.465 


8 








15.421 


-0.249 


8 


17.892 


-5.863 


7 








15.451 


-0.486 


8 


17.946 


-5.938 


6 








15.495 


-0.824 


7 














15.565 


-1.880 


7 














15.669 


-1.942 


2 















''Logaritlmi of the rest-frame frequency (Hz) in the middle of the bin. 

''Logarithm of the ratio of the monochromatic luminosity at the given frequency to that at 5000A. 
"^Number of photometric data points contributing to that frequency bin. 
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Fig. 1a. — SEDs for the first seven BAL quasars listed in Table 1; data are listed in Tables 2 and 3. The composite SEDs (dotted from 
Elvis et al. 1994, dot-dashed from Richards et al. 2006) have been normalized to the 24/im data point. The predicted X-ray SED (solid 
line) is normalized from the 2500A monochromatic luminosity (sec §4). The dashed curve is a model starburst SED normalized to the 
far-infrared photometry in excess of the composite quasar SED (see §4.3). Objects are labeled with name, BAL type, and redshift. 
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Fig. 1b. — See Figure la caption. 
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Fig. Ic. — See Figure la caption. 
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Fig. Id. — See Figure la caption. 
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Fig. 1e. — See Figure la caption. 
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Fig. If. — See Figure la caption. 
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Fig. 2. — Left: Histograms of cio values for the BAL (solid) and SDSS comparison (dashed) quasar samples. More negative values of 
ao indicate a redder ultraviolet-optical continuum. Right: A plot of ao vs. 'j^qq^ for the BAL (open diamonds) and SDSS comparison 

(filled circles) samples. Known LoBAL quasars are indicated with double diamonds. Both samples show a tendency for redder quasars to 
be more luminous at 5000A; this is likely a selection effect. 




Fig. 3. — In both panels, the BAL quasar (open diamonds; solid line) and SDSS comparison (filled circles; dot-dashed line) samples are 
plotted. Known LoBAL quasars are indicated with double diamonds. Best-fitting linear fits are drawn for each sample, and slopes are 
labeled. Left: Plot of Isiim vs. ^gggg^- The BAL quasar and SDSS comparison samples have both consistent slopes and normalizations 

(17.1±4.1 and 18.9±4.3, respectively) in this parameter space. Note that the la normalizations for both samples cannot be plotted within 
the bounds of this figure. Right: Plot of 'gQgQ^/'Sfim vs. 'g^Qp^i this representation gives a better idea of the relative luminosity of the 
optical versus the mid-infrared. Again, both slopes and normalizations are within la for the BAL and non-BAL quasar samples, and the 
la normalizations for both samples are outside of the figure bounds. 
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Fig. 4. — The far-infrared (1mm to 20fim) integrated luminosity in star formation, Lfir.sF vs. quasar far-infrared luminosity, iflr,qsO- 
Double diamonds indicate LoBAL quasars. Upper limits in the lower right of the plot were set to 20% of Lfi,. qsq when the data were 
consistent with the composite quasar SED. Note that the apparent gap in the distribution is from the disparate data available; all of the 
quasars with tight upper limits to ifir.SF have SCUBA data except for 2111—4335, which has a constraining 70/im upper limit. Presumably 
with complete submm coverage the distribution would more evenly fill this parameter space. 
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Fig. 5. — Mid-infrared-to-X-ray BAL quasar composite SEDs. The algorithm for constructing the composite SEDs is described in §4.4. 
The data and the Richards et al. (2006) SDSS composite SED (black curve) have been normalized to ^'^qqq^- ^'^^ reference, the 5000A 
normalization point (log{u) = 14.78) is marked with an arrow, and 8 fim is at log(!.') = 13.57. Top panel: All 38 BAL quasars (green 
filled circles) and their composite (black crosses); the SDSS comparison sample composite (blue diamonds) and its corresponding Richards 
et al. (2006) SED (dotted curve) are offset vertically for clarity. Bottom panel: The HiBAL quasar composite SED (black crosses) and 
the BAL quasar data delineated by BAL type: Hi (cyan x), unknown (green filled circles), and Lo (red *). 



